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r^ Abstract. In this paper we shall discuss the effect of disorder induced configuration fluctuations 

on single particle and two-particle phonon Green functions in substitutional random binary alloys. 
The randomness of the system will be dealt within the augmented space theorem, introduced by 
one of us [1]. This will be combined with a generalized Edwards-Langer diagrammatic technique 
to extract various useful results in the form of mathematical expressions. The purpose of present 
manuscript will be to describe the scattering diagrams in full detail. This will form the basis of the 
calculational algorithms which we had used in an earlier paper on NiPd and NiPt alloys [2, 3]. We 
shall show the structure of all possible scattering diagrams up to the 4-th order and subsequently 
illustrate how to obtain Dyson's equation from a re-summation of the diagrammatic series. We 
shall also study how disorder scattering affects two-particle Green functions associated with thermal 
response. We shall show that disorder scattering renormalizes both the phonon propagators as well 
as the heat currents. These renormalized heat currents will be shown to be related to the self-energy 



w _ of the propagators. We shall also study a different class of scattering diagrams which are not related 



to the self-energy but rather to the vertex corrections. 
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£T) ' PACS numbers: 72.15.Eb, 66.30.Xj, 63.50. 
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1. Introduction 

o 
\o . 

In an earlier communication [2] we had discussed the calculation of the configuration averaged thermal 
conductivity and diffusivities of disordered NiPt and NiPd alloys. We had argued that lattice 
thermal conductivity yields valuable information about the interactions of thermal excitations with 
composition fluctuations on the crystal lattice. We had also argued that although the theory of lattice 
thermal conductivity for perfect crystals and ordered compounds has been developed, the same is not 
true for disordered alloys. The presence of disorder results in scattering that not only depends on the 
impurity concentration but also crucially on both the relative masses and size difference between the 
O ■ constituent atoms as well as the difference between their dynamical matrices. For large mass and/or 

dynamical matrix differences, the effect of disorder can be quite unusual. That earlier communication 
focused on the results of numerical calculations on NiPt and NiPd and their comparison with available 
experimental data. In this communication we shall present a detailed analysis of the effects of disorder 
scattering on the thermal conductivity. We shall use the Edwards-Langer scattering formalism [4, 5] 
coupled with the augmented space technique introduced by our group [6] to analyze the disorder 
effect. This analysis will form the basis for justification of the calculational algorithms used by us in 
our earlier paper. 

X corresponding author, email : alam@bose.res.in 
f email : abhijit@bose.res.in 
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2. Effect of disorder on the single particle Green function : the Dyson equation 

The augmented space formalism (ASF) for carrying out configuration averaging of physical properties 
of disordered systems has been described in detail in several earlier papers (see [1],[3],[6]). We shall, for 
the sake of completeness, describe only those features which will be necessary for the implementation 
of our ideas in this communication. 

Let f(n R ) be a function of a binary random variable n R , whose binary probability density is given 
by: 

Pr(n R ) = x A S(n R ) + x B S(n R - 1) 

Such random variables are useful for describing randomly substitutional binary alloys. These are the 
kind of disordered systems we shall be interested in for this communication. The ASF now prescribes 
that we write this probability density, which is a positive definite integrable function as the resolvent 
of an operator whose spectrum consists of the random values taken by it (in this case and 1) : 



Pr(n R ) = -- %m(U\(n R I - Nr^^U) 

7T 



(1) 



The configuration space of n R is of rank two and spanned by the states |0) and |1). The operator N R 
acts on this space. | t/f) = \/xa\0) + \/xb\1) is called the reference state. Its orthogonal counterpart 
is I Ir) = \/xb\0) — \/x~a\1). The representation of N R in this new basis is : 

N _ ( %A y/x~AXB 



y/XAXB X B 

The ASF now proceeds as follows : 



</(///,.)»= / f(n R )Pr(n R )dn R = -i^m / /( /, „)•' | „ j i n „ I - X,. ) " '' | j ,<),!,, , 



i r°° 

= __ cj TO J2 J2 / /(^)(tfl|A)(A|(n J? /-7V fi )- 1 |A')(A'|t K )dn i? 

n A=0,1A'=0,1 °° 

= E (Tfl |A)/(A)(A| U) = (U |f| T«> (2) 

A=0,1 

Here f is an operator built out of f(n R ) by simply replacing the variable n R by the associated operator 
X R . The above expression shows that the average is obtained by taking the matrix element of this 
operator between the reference state | Ir). The full Augmented Space Theorem is a generalization of 
this for functions of many independent random variables {n R }. 

The theory of phonons consists of solving a secular equation of the form : 

(Mw 2 -D) u(R,w) = 

where u a (R, w) is the Fourier transform of u a (R, i), the displacement of an atom from its equilibrium 
position R on the lattice, in the direction a at time t. M is the mass operator, diagonal in real-space 
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and D is the dynamical matrix operator whose tight-binding representations are : 

M = E mR Sa P Pr 
r 

D = e( E <r] ft+E E *& *w 

ft [ R'jtR ) R R'^R 

(3) 

where the sum rule has been incorporated in the first term of the equation involving D. 

Here P R is the projection operator |-R)(i?| and T RR i is the transfer operator \R)(R'\ in the real space 

Ti spanned by the tight-binding basis {|-R)}- R,R' specify the lattice sites and a,f3 the Cartesian 

directions. m R is the mass of an atom occupying the position R and & RR i is the force constant 

tensor. 

We shall be interested in calculating the displacement-displacement Green matrix G(R, R' , w 2 ) : 
G(-R, R',w 2 ) = (R\ (Mio 2 - D) _1 \R') 

Let us now consider a binary alloy A x B y consisting of two kinds of atoms A and B of masses ttia and 
tub randomly occupying each lattice sites. We wish to calculate the configuration-averaged Green 
matrix <C G(R, R' ,w 2 ) ^>. We shall use the augmented space formalism to do so indicating the main 
operational results here. For further details we refer the reader to the above monograph [6] . The first 
operation is to represent the random parts of the secular equation in terms of a random set of local 
variables {n R } which are 1 if the site R is occupied by an A atom and if it is occupied by B. 

In terms of these, the mass operator can be written as : 

M = E m B + ur {5m) Sap Pr ; dm = mA - ms (4) 

R 

According to the augmented space theorem, in order to obtain the configuration-average we simply 
replace the random variables ur by the corresponding operators A^ associated with its probability 
density given by equation (1), and take the matrix element of the resulting operator between the 
reference states. For a full mathematical proof the reader is referred to [6] . 

n R — > N R = x I + (y- x) p R + y/xy T^ 1 

Similarly the random off-diagonal force constants $ RR i between the sites R and R' can be written 
as : 



K% = ® a AA n R n R> + $ bb(! " ««)(! - »«') + ® a AB n R {\ -n R >) + n w {l - n R ) 



^BB + \^AA + ^BB Z ^AB 



n R n R < + ($ab-$bbJ (n R + n w ) 

(5) 



The augmented space theorem [1] yields the configuration averaged Green function[3] as 

« G(R, R', w 2 ) »= ({0} ® R\ (g- 1 - Di) _1 1{0} ® R'), (6) 
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where the VCA green matrix g is given by 

g=(<M>cj 2 -«D>)-\ (7) 

with the Mass operator M and the Dynamical matrix operator D in the augmented space has the 
form 

M = A(m) I® I + B(m) J2 p R ® Pr + F ( m ) H T R ® Pr ' 

R R 

= < M > + M', 
where 

A(X) =« X » = (xX A + yX B ), 
B(X) = (y-x) (X A - X fl ), 
F(X) = ^xy(X A -X B ). 

x,y being the concentration of the components A and B of the alloy A x B y . p R and T R arc the 
projection and transfer operators in the configuration space describing the statistical behaviour of the 
system. 

< d >= - y^ J2 <c ® rr ' ^t^pr+YI XI <c ® rri > ^® T ««' ( 8 ) 

R R'jtR R R'yiR 

Pji and Trri are the projection and transfer operators in the real Hilbert space. 



< $,Rfl' >= x 2 $ AA + y 2 <5>bb + 2xy <f> AB . 
$'s in the right hand side are the force constant tensor between different combinations of atoms. Also 



D i = J2 { - Tr - 2 * rr ' \ ® p fl+EE ***- ® Trr> 

R [ R'=£R. J R R'^R 



with 



2 U , T?(m\ ,.,2 n-U 



T R = B(m) u 2 p R + F(m) ^ 7^ 






where 



D (1) = (»-it)*(i) 
D< 3 > = (j, - x) 2 $ (2) 

D (5) = xy $ (2) . 



c# 


+ T U ) + D (3) 

+ J ,R' J +1J RR' 


PrVri 


+ D 


(5) r U -7-T1 
RR' 1 R 1 R' ' 




D (2) : 


= Vxy $(i), 




D (4) 


= V^y (y - x) 
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Figure 1. The scattering vertices for the averaged Green function 



and 



$ (1) = x $ AA - y <5>bb + (y~ x)§ AB , 



$ 



(2) 



®AA + <5>BB - 2®AB- 



Physical interpretation of the vertices is transparent if we consider the | |fl) state to be the ground 
state or a state with no configuration fluctuations and | In) to be the state with one configuration 
fluctuation at the site R. The operator p R then measures the number of fluctuations (0 or 1) and 
Tp either creates or destroys a fluctuation at the site R. From equation (8) then we can interpret the 
following : 

The scattering vertices B,F and D^'-D*- 5 ' can be represented as shown in Figure f . The vertex F 
is diagonal in real space and causes a configuration fluctuation at a site due to disorder in the mass. 
Similarly, the vertex B, also diagonal in real space, counts the number of fluctuations at a given 
site. The remaining vertices labelled D are all off-diagonal in real space. The vertex D^ counts 
the number of fluctuations at one or the other of the sites associated with the vertex. D^ 2 ) causes a 
configuration fluctuation, again at one of the sites associated with the vertex. D*- 3 -* counts the number 
of fluctuations at both the sites associated with the vertex. D^ 4 ) counts the number of fluctuations 
at one site and causes a fluctuation at the other. Finally, D^ 5 ' causes a fluctuation at both the sites 
associated with the vertex. 



Equation (6) can also be expressed as 

« G(R, R', w 2 ) »= ({0} R\ ( g + g C-! g + g Dx g B-! g 



|{0}®#> 



(10) 
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The first term in Equation (10) gives : 

m®R\z\m®R') = s(R,R\w 2 ) 

In scattering diagram formalism, wc shall associate a propagator represented by a horizontal arrow 
for each factor g. 

The second term in equation (10) yields zero since ({0} i?|Di|{0} R') = 0. 

The third term gives : 

({0}0i?|gf3 lg D lg |{0}0i?') = 

S 1 S 2 S 3 S i {C} {C} 

({C} S 2 \ g \{C'} S 3 )({C'} S 3 | D! |{0} S 4 ) ({0} S 4 \ g |{0} R') 
A little more algebra yields the following contribution : 

({0}0 J R|gD 1 gf3 1 g|{0}0i 1 ") = 

J2 s(R,S!,w 2 ) (Fw 2 )g(S 1 ,S2,w 2 )6(S 1 -S 2 ) (Fw 2 )g(S 2 ,R',w 2 ) + 
SiS 2 

Y, E s(R,S!,w 2 ) (Fw 2 ) g(S u S 3 ,w 2 ) 6(Si - S3) T>% 4 g(S 4 ,i?> 2 ) + 
Si s$ s& 

Y, E 8^1 V) ( F ^ 2 ) s(SuS 3 ,w 2 ) 5(S 4 - S 4 ) T>% 4 g(S 4 ,i?> 2 ) + 

Si S3 S4 

Y E Sto Si,w 2 ) &sk s(&, S 4 , w 2 ) <5(5x - S 4 ) (Fw 2 ) g(S 4 , R', w 2 ) + 

Si S2S4 

Y E siR.S^w 2 ) T>% 2 g(S 2 ,S 4 ,w 2 ) S(S 2 - S 4 ) (Fw 2 ) g(S 4 ,i?> 2 ) + 

Si S2S4 

]T ]T s(R,S u w 2 ) Dg^ g(S 2 ,S 3 ,w 2 ) <5(5x - S 3 ) Dgk g(S 4 ,i?> 2 ) + 

SlS-2 S3S4 

J2 E sO^V) D S& S(^,S 3 , ^ 2 )<5(S 1 - S 4 ) Dg^ g(S 4 ,i?> 2 ) + 

S1S2 S3S4 

J2 E 8^1 V) Dg S2 g(S 2 ,S 3 V) <5(S 2 - S 3 ) D^ g(S 4 ,i?> 2 ) + 

S1S2 S3S4 

E E s(R,Si,w 2 ) Dg Sa g(S 2 ,S 3 ,™ 2 ) <5(S 2 - S 4 ) Dg S4 g(S 4 ,i?> 2 ) + 

S1S2 S3S4 

J2 E 8(^, ^, w 2 ) Dg } 5 2 8(^2, S 3 , w 2 ) 5(S 2 - S 4 )*(5i - S 3 ) Dg S4 g(S 4 , i?', W 2 ) + 

5iS2 £3 £4 

2 53 g(i?,S 1)W 2 ) Bf iS2 g(S 2 ,S 3 ,w 2 ) 6(S! - S 4 )S(S 2 - S 3 ) Dg 54 g(S 4 ,i?> 2 ) 

S1S2 S3S4 

(11) 

Referring to the above equation, we shall now build up the scattering diagrams. We have already 
associated scattering vertices with the terms in Di. These are the seven different types of scattering 
vertices as shown in Fig. 1. The dashed lines with arrow are associated with the delta function. 
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Figure 2. The possible scattering diagrams for n=2. 



A connected diagram to order 'n' is then built up by stringing together (n+1) propagators associated 
with g, connected by n-vertices with all fluctuation lines combined in pairs. For n=2 there are eleven 
possible scattering diagrams, whose algebraic form is given in Equation. (11). These diagrams are 
shown in Fig. 2. For n=3, the possible kinds of scattering diagrams are shown in Fig. 3. Note 
that it involve terms with contribution from B, D^ 1 ', D*- 3 ^ and D^ 4 ' as well. These scattering vertices 
can not sit either in the leftmost or in the rightmost positions, because then one of the associated 
pscudo-fcrmion green function line vanishes. In Fig. 3(A-F), each of the sections represents different 
distinct classes of diagrams. We call each sections as the topologically distinct scattering diagrams. 
There exists three more sections in addition to these (i.e. A-F) which are obtained simply by applying 
the reflection operators vertically to the diagrams in sections B, D and E. Hence in total, there exists 
9 topologically distinct classes of scattering diagrams for n=3. 

For n=4, there are various classes of diagrams as shown in Fig. 4. In this figure we have not shown 
the scattering diagrams with all sorts of combinations of fluctuation lines, but label the multiplicity 
within the brackets () for each diagram which actually gives the number of possible ways of combining 
the fluctuation lines for that particular diagram. For the diagram in the last column of Fig. 4(A), the 
middle decoration indicates a collection of all those diagrams lying in between two vertices such that 
a single pscudo-fermion line goes out ( from the left vertex ) and a single pseudo-fermion line comes 
in ( at the right vertex ) for n=4. Similarly for the diagrams in the last column of Fig. 4(E), the 
middle decoration indicates a collection of all those diagrams lying in between two vertices such that 
either (i) one pseudo-fermion line goes out (left vertex) and two pseudo-fermion line comes in (right 
vertex) or (ii) two pscudo-fcrmion line goes out (left vertex) and one pseudo-fermion line comes in 
(right vertex). While for the diagrams in the last column of Fig. 4(1), the middle decoration stands 
for a collection of all those diagrams lying in between two vertices such that two pseudo-fermion line 
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Figure 3. The possible scattering diagrams for n=3. 
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Figure 4. The possible scattering diagrams for n=4. 
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goes out (from the left vertex) and two pseudo-fermion line comes in (at the right vertex). Each of 
these diagrams in the last column are said to be the topologically distinct classes of diagrams. Hence, 
overall there exists nine topologically distinct classes of diagrams for n=4 also. 

Amongst the above sets of diagrams ( for n=3 and n=4 ), there exists three basic subsets namely 

• A set of separable diagrams, which are those that can be broken into two along a propagator line 
without also breaking a pscudo-fcrmion line. 

• A set of non-separable, non-skeleton diagram which can not be broken into two, but the VCA 
green function in this diagram can be renormalized. 

• Skeleton diagrams involving all the crossed and complicated structured diagrams. 

If we club together the contribution of all the skeleton diagrams calling this the self-energy, and allow 
all phonon Green functions except the left-most to be renormalized by the separable and non-separable, 
non-skeleton diagrams, we get the Dyson equation : 

<G>= g + gS«G» 

For homogeneous disorder we have shown earlier that we have translational symmetry in the full 
augmented space [7] . We can then take Fourier transform of the above equation to get : 

«G(q,£)» = g(q,£)+g(q,£)£(q,£) «G(q,£)» 




Figure 5. Topological structure of the skeleton diagrams for the self-energy S. The central dark 
semicircle represents all possible arrangements of scattering vertices to all orders. 
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The diagrams for the self-energy are skeleton diagrams all of which have the structure as shown in 
Fig. 5. Each of these diagram starts or ends with any one of either F, D^ 2 ^ or D^ 5 ) vertex. The 
central dark semicircle stands for all possible arrangements of scattering vertices to all orders. 

3. Effect of disorder on two-particle Green functions : Lattice thermal conductivity 

3.1. Configuration averaging of lattice thermal conductivity 

The Kubo formula which relates the optical conductivity to a current-current correlation function is 
well established. The Hamiltonian contains a term ^jj« ' A(r, t) which drives the electrical current. 
For thermal conductivity we do not have a similar term in the Hamiltonian which drives a heat current. 
The derivation of a Kubo formula in this situation requires an additional statistical hypothesis which 
states that a system in steady state has a space dependent local temperature T(r) = [ftB/3(r)] _1 . The 
expression for the heat current has been discussed in great detail by Hardy [8] and Allen and Feldman 
[9] . The readers are refereed to these papers for the details of calculation. The matrix element of the 
heat current in the basis of the cigenfunctions of the Hamiltonian is given by : 

S^(k) = \ (u, k7 + u; ky ) < 7 ,(k), (12) 

where, the phonon group velocity v 77 / (k) is given by 



'--VS=*7??« t '« (k >l£^ H""" 



2 A /^k7^k7 



l== ££<(k) V k £"-(k) e ;,(k), 



fl V 



(13) 



here i,j label the atoms, 7, 7' label the various modes of vibration, fi, v label the Cartesian directions, 
Wk 7 ,Wk 7 ' are the eigen frequencies, e^(k),ey(k) are the polarization vectors associated with the 7-th 
and 7'-th mode of vibration. D M!/ (k) is the Fourier transform of mass scaled dynamical matrix. 

We shall consider the case where the temperature gradient is uniform within the system. The linear 
heat current response is related to the temperature gradient field via a generalized susceptibility : 

/CO 
dt! n^it-t') V v 5T{t), 
-°° 

The generalized susceptibility is the thermal conductivity in this case and 
«"» - 6(r) I / d\(S»(-ih\),S»(T)), 



Tjo 



0(t) is the Heaviside step function, and 
S(-ih\) = e XH S e- XH . 
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( ) on the right-hand side stands for thermal averaging over states in the absence of the temperature 
gradient. The above equation can be rewritten in the form of a Kubo-Greenwood expression 



•n«.T>-£/£ £ 



rfj(u),T) 



k b T 2 



7,7'#7 

8tt 3 



|g™s^( k) s;, 7 ( k )^ k7 -. ky - w ) 



(14) 



]T(n k7 ) S^k) 



rf 3 k 
8^ 



^(n k7 ) S^ 7 (k) 



d 3 k- 



- k W^E 



<9(n k7 ) 



87r 3 "^ <9(?ia> k7 ) 



5^ 7 (k) 5^(k) 



*(w) 



(15) 



where (n k7 ) = (e"™ u ' k 7 — 1) _1 is the equilibrium Bose Einstein distribution function and T is the 
absolute temperature. For an isotropic response, equation. (14) for inter-band transition can be 
expressed as 

«/(«, T) = Jl J^ J **' J £ § 77' ( k > T ) § 7'7( k > T )<V - ^ )*("' + w - c k7 ) 



where 



S^(k,T) 



(nky) - (n k7 ) 



We may rewrite the above equation as 



S^y(k). 



Ki(u>,T) 



' £ W^T 



3ttT 



8tt 3 



S"(k,T)Sm{G(k,w')} S"(k,T) 3?m{G(k, u/ + w)} 



(16) 



The operator G(w) is the phonon Green operator (Mu> 2 I — 3?) _1 . The Trace is invariant in different 
representations. For crystalline systems, usually the Bloch basis {|k, 7)} is used. For disordered 
systems, prior to configuration averaging, it is more convenient to use the basis {|k, a)}, where k is 
the reciprocal vector and a represents the coordinate axes directions. We can transform from the mode 
basis to the coordinate basis by using the transformation matrices T 7Q (k) = e 7 (k). For example 



S^(k,T) = T- 7 T (k)S^ 7 ,(k,T)T-yk). 



If we define 



k(z 1 ,z 2 ) 



fk 

8tt 3 



Tr 



SG(Mi) S G(k,z 2 ) 



(17) 
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13 



12ttT 



k"" (u/ + , u/ + + w) - k/" 1 (w'- , J' + lo) 



(18) 



where 



/(w*) = lim /(w ± W). 



We have also used the herglotz analytic property of the Green operator 

G{v + i8) = $te[G{u>)] - % sgn(<5) 9m [G(w)] . 

For disordered materials, we shall be interested in obtaining the configuration averaged response 
functions. This will require the configuration averaging of quantities like k(^ 1; z 2 )- For disordered 
materials, Equation. (18) should be expressed as 

<kj(w,T)» = — !— V /"dw' ((\k^(w'-,u;' + + w) + K> i > i {uj' + ,u>'- +w) 

12irT z — ' / \ \ L 

-«/"V + ,u>' + + u) - k""(u'" ,«'- + w)] )), (19) 

Let us discuss the configuration averaging of the two particle green function of the kind k(zi, zi). The 
augmented space theorem immediately implies that 



« k( Zi ,z 2 ) »= Tr ({0} SG(zi) S^G(z 2 ) {0} 



(20) 



The first thing to note about Equation. (20) is that the right hand side is an average of four random 
functions whose fluctuations are correlated. The average of the product then involves the product of 
the averages and other contributions which come from averages taken in pairs, triplets and all four 
random functions. 

Following a similar procedure as for a single particle green functions, the operator S in the augmented 
space has the form 



-11 



S = EE[« § »Ro,w a ' I + Sj£ ffl , a , (p R + Pr ) + S^lnw fa 1 + ?l 

Ra R'a' 

^rIr'^Pr^Pr* + 8%L,r.*> (pr ® T$ + p R , T^) + S%i Rla ,T^ T# 



where 



®T RR ,, 



(21) 






xyS {1) , 
xy (y-x) S (2) , 
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Figure 6. The scattering vertices associated with the random current terms. 
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3.2. The disorder renormalized current 



We now start to set up the scattering diagrams for the thermal conductivity. A look at Equation. (21) 
shows us that the first term is the averaged VCA current. This term is absorbed in the unscattered 
part of the phonon green function and leads to the zero-th order approximation. Equation. (21) looks 
very similar to Equation. (9) from the operator point of view. The only difference is that the former 
equation arises due to the disorder in heat currents while the latter equation due to the disorder in 
the dynamical matrix. Exactly as before we can associate scattering vertices with the terms in S. 
Fig. 6 shows sixteen different scattering vertices arising out of the equation (21). Let us now discuss 
how the scattering diagrams are set up and then examine them. The rule for obtaining the diagrams 
for the correlation function <C k(zi, zi) 3> is as follows : Take any two current diagrams from Fig. 6 
and two propagators and join them end to end. Now join the configuration fluctuation lines (shown 
as dashed arrows) in all possible ways. The zero-th order approximation for <^ n(zi,Z2) 3> can be 
shown diagrammatically as in Fig. 7. The most dominant contribution comes from this particular 
diagram. Here the two current terms are the averaged current, and all configuration fluctuation 
decorations renormalize only the two phonon propagators. The bold propagators in this diagram are 
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Figure 7. The VCA or zero-th order approximation for <C k(zi,Z2) 2>- 

fully scattering rcnormalizcd propagators corresponding to the configuration averaged green function. 
The contribution of this term to the correlation function <C k(zi, z-i) 3> is 



8^ 



« S(k) X G(k, Zl ) X S(k) >< G(k, z 2 ) > 



(22) 



The rest of the terms in Equation. (21) give rise to scattering. We shall now focus on the main 
correction terms to the expression in Equation. (22). These are the correction terms to the averaged 
current which, as we will show, are closely related to the self-energies. The first type of scattering 
diagrams are those in which no disorder propagator (shown as the dashed lines) joins either two 




Figure 8. Few examples of scattering diagrams where no disorder line joins the two phonon 
propagators. 
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phonon propagators or two of the current lines directly. Fig. (8) shows few such scattering diagrams. 
These sets of diagrams may be clubbed together and renormalized in a form which will consist of two 
fully renormalized phonon propagators connected at the two ends by a new form of the renormalized 
current. This new form of the renormalized current may be obtained in the following way. 

Fig. 8 clearly shows that these types of diagrams are made out of a left renormalized current diagram 
chosen out of any one of the diagrams from (a,b) in Fig. 9 and one right renormalized current diagram 
from any one of the diagram (c,d) of Fig. 9 connected by two renormalized propagators. 




(C) (d) 

Figure 9. Scattering diagrams contributing to effective heat current. 
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Let us now obtain expressions for the renormalized currents. A careful look at the self-energy diagrams 
(See Fig. 5) shows that all self-energy diagrams have the structure 



£(k,*) = A(k,z) J7(k,z) A(k,z), 
Where $~2(k, z) is the Fourier transform of 

JW(*) = E G ««^) P h!rM Gr 2 r>{z), 

-R1-R2 



(23) 



(24) 



and 



A(k, z) = F z 2 + 2 D (2) (k) + 2 D (5) (k) 



In the above equation, the quantity 'P' stands for the central dark semicircle of Fig. 5 which represents 
all possible arrangements of scattering vertices to all orders. 

If we compare the diagrams of Fig. 9(a) with the diagrams for the self energy Fig. 5, we note that 
the only difference between the two is that the left most scattering vertex is replaced by a very similar 
current term. In the diagrams of Fig. 9(a), the left most diagonal terms similar to the vertex F of 
Fig. 1 is off course missing. The contribution of such diagrams may be written in a mathematical 
form as 

(2 S( 2 '(k)+2S (5) (k)) n(k,z) A(k,z). 
which may be expressed in terms of the self energy 'X' by using Equation. (23) as 



£(Mi). 



(2S( 2 '(k) + 2S< 5 )(k)) [A(k, Zl ) 
The contribution of the diagrams in Fig. 9(b) is 

£(k, z 2 ) [A(k, z 2 )] ~* (2 S( 2 ) (k) + 2 S^ (k)^ 
Similarly the contribution of the diagrams in Fig. 9(c) and Fig. 9(d) are respectively given by 

S(k, Zl ) [A(k,z 1 )]' 1 (2S( 2 )(k)+2S( 5 '(k)) & (2S< 2 )(k) + 2S< 5 )(k)) [a(M 2 )] -1 S(k,z 2 ) 

The next most dominant disorder corrections come from a group of diagrams which describe joint 
fluctuation of one current and two propagators. A few such diagrams are shown in Fig. 10. 

The contributions of these diagrams [Fig. 10(a,b)] can also be expressed in terms of the self energy 



as. 



(a) S(k,z 2 ) fA(k,« 2 )l S< 5 )(k) [A(k,«i)l E(k,*i), 



(b) S(k, Zl ) [A(k,«i)l S^(k) [A(k, 



Z2J 



S(k,z 2 ) 



If we now gather all the contributions from these diagrams [ From Fig. 9(a,b,c,d) and Fig. 10(a,b) ], 
we may define a renormalized current term as follows : 



S cff (k,zi,;z 2 ) =< S(k) > +ASi(k,zi,z 2 ) + AS 2 (k,zi,z 2 ), 



(25) 
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Figure 10. The scattering diagrams associated with joint fluctuations of one current term and two 
propagators. 



where 



ASi(k,z!,z 2 ) = 2 (S^(k) + S^(k)j |A(k,zi)J S(k,«i) 
S(k, z 2 ) f A(k, z 2 )l ~* 2 fs^(k) + S( 5 )(k)) , 



AS 2 (k,z 1 ,z 2 ) = E(k,z 2 )[A(k,z 2 )] S< B >(k)[A(k,*i)] E(k,*i). 

The contribution of these disorder-renormalizcd currents and propagators to the correlation function 
is 



/d 3 k r 
^-3 Tr |S eff (k, zi,z 2 ) < G(k, z x ) > S^ ff (k, zi,z 2 ) < G(k, z 2 ) > 



(26) 



We shall now discuss the disorder correction terms which involve joint fluctuations between the two 
current terms and one propagator. Few such diagrams are shown in Fig. 11. A close inspection of 
these diagrams shows that these are also related to the self energy diagrams with vertices at both ends 
replaced by currents. The corrections due to these terms can therefore be related to the self-energy 
as before. 



The contribution of these diagrams to the correlation function is given by 
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Figure 11. The scattering diagrams associated with joint fluctuations of two current terms and one 
propagator. 



where 



/d 3 k 
^ Tr [AS 3 (k, zi) < G(k, zi) > + AS 4 (k, z 2 ) < G(k, zi) >] 



(27) 



AS 3 (k,«i) = (sW(k) + S^(k))[A(k,«i)] E(k,«i) [A(k,«i)] (s< 2 >(k) + S^(k)J , 

AS 4 (k, z 2 ) = (S( 2 ) (k) + S< 5 )(k)) f A(k, z 2 )l ~* E(k, z 2 ) f A(k, z 2 )l ""' f S^ 2 ) (k) + S< 5 )(k " 



In our earlier paper [10] on a similar problem, we have argued that these are the dominant disorder 
corrections to the average current. Intuitively we also expect the same to be true in the present case 
as well. It is important to note that that these corrections can be obtained from the self-energy and 
is therefore eminently computationally feasible in the case of realistic alloys, once we have a feasible 
method for obtaining the self-energy. 

There are other scattering diagrams which are not related to the self-energy but rather to the vertex 
corrections. In these diagrams, a disorder line connects both the phonon propagators directly We 
expect the corrections from these types of diagrams to be less dominant. For the sake of completeness, 
we shall indicate in detail how to obtain them within a ladder diagram approximation in the next 
section. 



Phonons in disordered alloys : a multiple scattering approach 



20 



Scattering diagram 



Vertex correction 





Maximally crossed 



fe, 



Effective 
Current 



Vertex 



Figure 12. The scattering diagrams leading to vertex correction. 

3.3. The vertex correction 

The vertex corrections are basically those scattering diagrams in which disorder lines connect both 
the propagators directly. We have not yet incorporated these kinds of diagrams in the disorder 
rcnormalization. These types of diagrams arise due to the correlated propagation. The diagrams 
leading to the vertex corrections may be of different kinds e.g. ladder diagrams, maximally crossed 
diagrams etc. The ladder diagrams are those diagrams which are built out of repeated vertices shown 
on the first line of Fig. 12. These kinds of diagrams can be summed up to all orders. This is 
the disorder scattering version of the random-phase approximation (RPA) for the phonon-phonon 
scattering. The maximally crossed diagrams are those diagrams in which the ladder inserts between 
the crossed vertices. These types of diagrams are shown in the second line of Fig. 12. Here we shall 
consider the ladder diagrams in detail and show how to obtain the contribution of these diagrams in 
terms of mathematical expression. We then sum these ladder diagrams to all orders. The possible 
scattering diagrams for the ladder kind of vertex correction involving the vertices B, F, D^ 1 ' to D^ 5 ' 
are as shown in Fig. 13. 

The contribution of seven categories ( A-G ), shown in Fig. 13, of the ladder scattering diagrams in 
terms of mathematical expression are as given below. 

Category A 



{ W $) A = (^) 2 F a F- t 5 a[i 5 1& + 2 4 D ( H F"> S jS 
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Figure 13. The ladder scattering diagrams for the vertex correction. 
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2 «) 2 



/ y {Fqv'^av') Grv',Rv" {F v » p5 v i 



a) 



+ (z 1 2 B afj 5 af3 )(z 2 2 B„ s 5 u5 ) + 2 \d { ^(z 2 2 B uS 5 u5 ) 



x{z 2 2 B vS 5 vS ) + 2 



2_^( Zl Fav'$av')GRv',Rv"D v ii 



(2) 



2 y D av , Gr u ^R u " (Zl FjylipS^n 



(B v $5 vS ) 

E D uv>Gris',Ru"D v , ii3 
v'v" 

{z 2 2 B v5 5 v5 ) 

(5)^2; 



{z2 t B vS 6 v s) + D^{z 2 2 B vS 5 vS ) 



Category C 



(w$) c = 2[( Zl 2 B a0 S af) )D% 



/ X Z ^ Fav'&av') Gr v '.b,ij" {z\ F v " fib~ v " p) 



D 



(1) 



4 D (1) D (1) 



Gru'^Ruii D v „p 



E ^ g*,,^ »$ 



(2) 



_v' v" 



V® +2 



D 



(i) 
i/5 



y^(zi 2 F a ^(5 a ^) 



_v'v" 



E D 3 Gr v ',r v » {zi 2 F u „p5 v „ 



D 



(i) 

v5 



Category D 



(w$) D = 2[( Zl 2 F af3 6 af3 )D% 



+ AD {2) D {2) 



Category E 



Ki) j 



= 2 



D 



(3) 



E^ 1 B av i5 av i) Gb.v',Rv" [Z\ B v ii p5 v 'i fi) 



2_^{ z i B au iS au i) 



Grv'^rv" D v „p 



D% +16 



D 



_v'v" 

(3) , 
i/5 "T" 



/ : 1^1 -^aa/i </cki/i J ^ Rv\.Rvi \Z\ " v%v 3 Qviv$ ) ^Rv^^Rvk \^1 ^1/41/5^1/41/5/ 
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Figure 14. The structure of infinite series of ladder diagrams contributing to the correlation function 
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Therefore, the sum of all possible scattering diagrams contributing to the four legged vertex ( shown 
in the extreme right column of Fig. 13 ) will be given by 



w$ = E K 7 






Here we shall sum the ladder diagrams to all orders. The contribution of a single ladder diagram to 
the correlation function as shown in the top line of Fig. 14 is 

/ j / , / , / j / j / j ^Rr a & M iai GR lCtu R 2a2 (z 1 ) W"^s & R 2 a 2 ,R 3 a 3 ( z l) 

-R1-R2 R3R4 -R5 ai«2«3«4CI5a6 

\^R 3 a 3 ,R4.aJ G Ria4 M 2aa (z 2 ) G R2 a 5yRsa6 (z 2 ) (28) 

If we apply the homogeneity in full augmented space, it will imply that the above expression is 
independent of 'R' which allows us to take the Fourier transform leading to 
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where we have defined 



T( Zl ,z 2 ) = / - T G(k,z 2 )S eff (k,^ 1 ,z 2 )G(k,z 1 ), 
J biz 8tH 

T(z u z 2 )= / -^G(k',z 1 )(S off (k',z 1 ,z 2 )) t G(k',z 2 ). 

Let us now look at the contribution of the infinite series of ladder diagrams ( shown in the third 
column of Fig. 14 ) to the correlation function. Each one of them has the same structure as Equation. 
(28), We may then sum up the series as follows : 

Let us define 

d 3 k 



®l S a(zi,z 2 ) = / — -3 G Q/3 (k,zi) G 7d -(k, z 2 ) 

J BIZ 87r 



IBIZ 

Then 

A(zi,z 2 ) = W + W0W + W0W0W + ... = W(zi,z 2 )(l-®(zi,z 2 )W(zi,z 2 j) . 

Thus the contribution of the infinite series of ladder diagram vertex corrections to the correlation 
function may be expressed as 

« AK(z u z 2 ? adder » = ^^r^(z 1; z 2 ) A%(Zl,Z2) T V S (Z1,Z2) 

3 fS 

\T{z u z 2 ) f (zi, z 2 ) A(«i, z 2 )] . (29) 



a/3 7(5 

Tr 



Including all kinds of disorder corrections, the configuration average of the two-particle Green function 
( or the correlation function ) now has the form, 

< k(zi,z 2 ) >=< k (1) (zi,z 2 ) > + < K(2)(zi,Z2) > + < AK(zi,z 2 ) lQdder > . (30) 

4. Conclusions. 

We have shown that the effect of disorder scattering is to renormalize the phonon Green functions, via 
a Dyson equation. Disorder scattering also renormalizes the heat currents. These corrections belong 
to two classes : one which is related to the self-energy and another related to the vertex corrections. 
Once these relations are established, we can use the augmented space recursion algorithm [3] to obtain 
the self-energy and the current corrections from it. The expression for the current corrections related 
to the vertex corrections has also been explicitly shown here. These will form the justification for the 
calculational algorithm which we had used in our earlier paper on NiPd and NiPt alloys [2]. 
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